Domain Formation: Decided Before, or After the Transition? 
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There is increasing evidence that causality provides useful bounds in determining the domain 
structure after a continuous transition. In devising their scaling laws for domain size after such a 
transition, Zurek and Kibble presented arguments in which causality is important both before and 
after the time at which the transition begins to be implemented. Using numerical simulations of 
kinks in 1+1 dimensions, we explain how the domain structure is determined exclusively by what 
happens after the transition, even though the correlation length freezes in before the transition. 
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I. INTRODUCTION 

It is inevitable that causality will provide constraints 
on the formation of domains after a continuous tran- 
sition. Although, adiabatically, correlation lengths di- 
verge, the finite speed at which the order parameter can 
adjust to the changing environment guarantees that cor- 
relation lengths remain finite for transitions implemented 
in finite times. Assuming that systems do change as fast 
as they can, Kibble [| and Zurek H Q predicted simple 
scaling behaviour for the maximum correlation length of 
the order parameters of the early universe and condensed 
matter systems respectively, as a function of the cooling 
rate. 

Domain formation, the frustration of the order param- 
eter fields, is often visible through topological defects, 
typically vortices, which mediate between different equiv- 
alent ground states. Since defects are, in principle, ob- 
servable, they provide an excellent experimental tool for 
confirming this scaling behaviour, insofar as the defect 
separation can be correlated directly to the order pa- 
rameter. Several experiments have been performed to 
measure defect density that support the KZ scaling be- 
haviour under this assumption 0. 0. Q ■ which we term 
the Kibble-Zurek scenario, or are commensurate with its 
predictions 0, 0,1111113 • 

In this paper we give a critical reappraisal of the 
Kibble-Zurek (KZ) scenario or, more accurately, scenar- 
ios. There are several different ways to derive causal 
bounds on defect densities. Although they agree nu- 
merically for simple systems, they differ conceptually in 
their assumption as to how the length scale that deter- 
mines the defect separation arises. The main distinction 
is whether this is set when the system freezes in before the 
transition or when the system relaxes after the transition 
has begun. 

In this paper we shall show that despite the fact that 
the correlation length freezes in before the transition, this 



value does not determine the density of defects upon 
their appearance, contrary to the standard formulation 
of causal bounds. As we shall see, the defect density is 
determined entirely by what happens after the transition 
has begun. This latter point has been argued by one of 
us (RR) over many years 0, 0, 0, 0] in the context 
of analytic approximations and by Moro and Lythe [l5| 
using numerical simulations to supplement similar ana- 
lytic results but, we feel, can benefit from the explicit 
demonstration that we give here. It is also timely in that 
the first experiments [16| are being performed that show 
that the defect density depends on what happens after 
the transition. We shall turn to these in our concluding 
section. 



II. THE KZ SCENARIOS 

To be specific, consider a system with critical tempera- 
ture T c , cooled through that temperature so that, if T(t) 
is the temperature at time t, then T(0) = T c . 

T(0) = ~T c /t q 

defines the quench time tq. The adiabatic correlation 
length £ a d,(t) = £ad(T(t)) diverges near t — but the 
true correlation length £(t) remains finite. The standard 
formulation of the KZ scenario is that the correlation 
length £ that characterizes the onset of order is the equi- 
librium correlation length £ = £ a d{i) at some appropriate 
time t < before the transition is implemented. 

The argument goes as follows; at time t there is a max- 
imum speed c(T(t)) = c(t) at which the system can order 
itself. In the early universe this is the constant speed of 
light, in superfluid 4 He, say, the speed of second sound, 
vanishing at T — T c . [Whether there is critical slowing 
down or not is irrelevant.] ^From this we can define a 
relaxation time r(t) = £ a d(t)/c(t). The earliest sugges- 
tion, due to Zurek, is that, initially, for t < some way 
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before the transition, the system behaves adiabatically. 
However, the system freezes in (the 'impulse' region) at 
time — i, where 

f~T(t), (1) 

i.e. when the time to the transition matches the relax- 
ation time. It is proposed that the correlation length 
£ = £ad{~t) frozen out at this time will determine the 
domain size at the transition. 

The argument due to Kibble is even simpler; can- 
not grow faster than c(t). If <^ (t) « £ a d(t) until the system 
cannot keep up with the cooling, the time —t at which 
the system freezes in is defined by the condition 

Ld(~t) » c(-i). (2) 

For simple systems both estimates of i agree, to give the 
allometric form Q 

U=U(-t) = a < (T Q y<. (3) 

The scaling exponent a < depends on the (adiabatic) criti- 
cal behaviour of the system as T — > T c — , and a< depends 
on the microscopic properties of the system. 

In addition to considering relaxation times before the 
event Zurek also invoked [3j causal horizons after the 
event to obtain causal bounds. At time t the causal hori- 
zon in which the order parameter can be correlated is 

&(i)«2 / dsc(s). (4) 
Jo 

^From this viewpoint our horizon bound for the earliest 
time i that we can see defects is when £ c (t) becomes larger 
than the coherence length £ a d{t). That is, the causal 
horizon is big enough to hold a classical defect, 

(5) 

The correlation length is then £> = £ a d(i). 

A parallel case that the relevant time for setting £ is 
after the transition can be made for Kibble's analysis, 
by saying that the impulse regime will give way to the 
adiabatic regime once 

Ld{i) « -c(t), (6) 

and it is this that determines the relevant t. Insofar as 
£,ad{t) is symmetric, this gives the same t as in J5J). Again, 
both estimates © and © for t agree, to give 

£">=a>(T Q r>, (7) 

where the scaling exponent er> depends now on the (adi- 
abatic) critical behaviour of the system as T — ► T c — . 

For continuous T(t) the scaling exponents cr< = cr> 
are equal, when both exist. As a result and J7J) are 
conflated as £ — £i(tq) ct , the KZ scaling law. 

In either formulation Kibble and Zurek then make the 
second assumption, that the relevant £ sets the initial 



scale of the defect network, with defect separation £def ~ 
I 

In this paper we shall show that despite the fact that 
the correlation length freezes in at £ satisfying ©, this 
does not fix the density of defects. Instead, the density is 
determined by the growth of instabilities after the tran- 
sition, which essentially matches jjj without invoking 
causality directly, although all evolution is, inevitably, 
causal. 

III. THE MODEL: LENGTH SCALES 

We will consider a simple model in 1+1 dimensions 
that has been used as a first step in several studies of 
defect formation [T7L Hsj . This model is easily amenable 
to numerical simulation and leads to the correct scaling 
laws for the defect density. In a future publication we will 
deal with systems with larger number of spatial dimen- 
sions. Nevertheless we expect the results discussed below 
and the analytical estimates behind these to generalize 
to higher dimensions. 

Specifically, we consider the following Langevin equa- 
tion, describing a time-dependent Ginzburg-Landau 
(TDGL) field theory: 

d 2 (j) - V 2 + a 2 d t 4> + m 2 (t)(j) - 2A</> 3 = a£, (8) 

where £ is a Gaussian noise term obeying 

<C(x', t')((x, t)) = 2TS(x' - x)5{t> - t), (C(x, t)) = 0. 

. (9) 

In one dimension V 2 — d 2 <p, a measures the amplitude 
of the noise, and its square quantifies the dissipative ef- 
fects. This relation between the amplitudes of the noise 
and the dissipative term ensure that the fluctuation- 
dissipation condition is satisfied and that for very long 
times the system should eventually reach thermal equi- 
librium at temperature T. 

This TDGL model is motivated empirically for over- 
damped (large a) condensed matter systems but, even 
in their idealised form, we would expect relativistic field 
theories to have non-linear noise arising from their envi- 
ronments. However, we have shown elsewhere that scal- 
ing behaviour is, essentially, independent of the linearity 
or non-linearity of the noise and we stick with JSJ. 

We start by presenting a simplified analytic argument 
supporting the assumption that the correlation length 
does freeze in before the transition, as anticipated by 
Kibble and Zurek. This will also help clarifying the rea- 
son why this scale does not necessarily set the value of 
the final defect density. In the simplest version of the 
model, the mass term is decreased linearly in time from 
an initial positive value /1 2 as 

m 2 (t) = -^ 2 — , -T Q <t<0 (10) 

For t < —tq the mass term is set to constant /z 2 . For 
values of time larger (smaller) than tq (— tq) the mass 
term is set to constant — /1 2 (/i 2 ). 
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To understand how the system does freeze in at time 
—i, for t < non-linear excitations have not had time to 
grow and it is sufficient to set A = in (0. For reasons 
of familiarity it is convenient to extend the (now linear) 
equation to three space dimensions for t < 0. The field 
correlation function G(r,t) = (<f>(x)(j)(0))t is expressible 



G(r,t) 



d 3 k 
(2^)3 



(ii) 



where P(k, t) is the power in the fluctuations of wave- 
vector k. The solutions to © are given in terms of Airy 
functions and their generalisations, and are messy with- 
out being very informative (see |15|). It is sufficient, for 
illustrative purposes, to restrict ourselves to the case of 
strong dissipation, whereby Airy functions become ex- 
ponentials. As we have shown before \v^, for strong 
damping, 

P(k, t)= dr T(t - r/2) e" Tfc2 e" Jo dT m (*~ T ^ , 
Jo 

up to an irrelevant renormalisation. In turn, this gives 
G (r,t) = 



oo / y \3/2 

drT{t-T/2) - — ) r • 
o V 47rr 



T/2)( — ) e - 2 /4r e -/>'^'/^ 

(13) 



It is not difficult to see that, for r large enough, 



G(r,0) 



1 

47rr 



exp[a(r/£<) 



4/3i 



(14) 



in this case. In (|14|l £< is, indeed, the correlation length 
that follows from Q and @ above and awl. 

What is important here is that the correlation length 
£< is determined from the position of the nearest sin- 
gularity in the complex fc-plane of the integral (|llfl after 
performing the angular integrations or, equivalently, from 
the large-distance behaviour of the correlation function. 

The reason why, both for this and the more general 
case of arbitrary a, this length does not set the scale for 
the density of defects is that defects are identified by the 
field zeroes at their cores when they are formed at t > 0. 
Initially, prior to t = 0, there are zeroes on all scales. 
However, once the transition begins to get under way, 
the instabilities in the long wavelength modes of the field 
lead to exponential growth in long wavelength amplitudes 
that orders the fields. Insofar that we can persist with 
the linearized system after the transition, the separation 
of zeroes t; zer o at the formation of defects is |2(], 0] 
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-l G"(r = 0,f) 
2^C(r = 0,(")' 



(15) 



where dashes denote differentiation with respect to r. 
That is, the information comes entirely from the short 
distance behaviour of G(r,t) for t > 0, rather than the 



long-distance behaviour that determined £< of (|14|) for 
t < 0. It is determined by the second moment of the 
power spectrum, rather than its analytic structure. 

It happens that, in mean- field theory at least, £ zero of 
(|15fl can have the same scaling behaviour as £> of J7J, 
that would be obtained from © and JSJ. There are only 
logarithmic corrections to the scaling laws, provided ther- 
mal fluctuations, that introduce a further length scale 
(the Ginzbur g co rrelation length) are unimportant [l3| . 
The work of [l5| shows just how reliable this analytic 
approximation is, provided transients are taken into ac- 
count. 



IV. NUMERICAL SIMULATIONS 

For this type of quench, we expect that some time after 
tq the field should settle to a series of alternating posi- 
tive and negative vacuum regions. These regions will be 
separated by the topological defects of the theory, kinks 
and anti-kinks interpolating between the opposite vacua. 
The defects will then enter a regime of slow evolution 
with pair annihilation taking place at long times. Dur- 
ing this stage of the evolution there is a one-to-one corre- 
spondence between kinks and anti-kinks and zeros of the 
scalar field <j), making it trivial to identify them numeri- 
cally. 

Depending on the value of the dissipation coefficient 
a 2 , we expect that the final density of defects N& c t = 
1 /£dcf w iU scale with different powers of the quench time 
scale. For high dissipation the system is effectively first- 
order in time, whereas for low a the second derivative 
dominates the evolution and the system behaves in a rel- 
ativistic fashion. Assuming that the final defect density 
is of the order of the inverse of the freeze-out correlation 
length, iVdof should scale as Tq° . As discussed exten- 
sively in the literature and verified numerically for this 
model [TtI ITsf , the exponent is given by a ~ 1/4 in the 
dissipative regime, and a ~ 1/3 in the relativistic case. 

We will now look at variations of the quench model 
described above in which we modify the basic system so 
that the properties of the quench are qualitatively differ- 
ent before and after the transition takes place at t = 0. 
By measuring the final defect densities and determining 
the corresponding scaling powers, we will be able to de- 
termine which period of the evolution dictates the out- 
come. Though transitions such as envisaged here are 
unlikely to be realisable in a physical situation, such sce- 
narios are useful as idealized experiments. We stress that 
our primary concern here is with a matter of principle, 
given the general acceptance of the KZ scenario in its 
original form. Nevertheless, as we will discuss below, ex- 
perimental results in transitions in Josephson Junctions 
can be understood in terms of distinct behaviours of the 
system before and after the transition. 

As already mentioned, the value of the dissipation in 
equation (jHJ determines the scaling power of the final 
density of defects. In the first variation we will rely on 
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FIG. 1: Final defect density scaling power as a function of 
the noise amplitude for t > 0. The amplitude before the 
transition was set to a< = 0.4 and q< = 0.5 for the full 
and dash-dotted curves respectively. For comparison we plot 
(dashed curve) the standard result for a case with constant 
dissipation throughout the quench, i.e. q< = q>. The error 
bars are the standard deviation of the result over a set of 10 
independent realizations of the quench. 



this fact to discriminate between the effects of the two 
main stages of the quench on the final outcome. The 
idea is to allow a to take two distinct values, a = a< for 
t < and a — a> for t > 0, characterizing the dissipa- 
tion and the thermal noise before and after the transition 
respectively. If a< and a> are values typical of differ- 
ent regimes (under and over-damped and vice-versa) for 
t < and t > respectively, the measured value of a 
should reflect the relevant period of the evolution. 

The numerical procedure used in this work follows very 
closely that of 01 • We will outline it here briefly and re- 
fer the reader to the aforementioned publication for a 
more detailed description. The parameters of the poten- 
tial are fj, 2 — 1.0 and A = 1.0 and the bath temperature 
T is set to a low value, typically T = 0.01. For every 
fixed choice of a< and a> a set of quenches is performed 
with tq = 2™, n = 1, 2, .., 9 and the final defect density 
is obtained by counting zeros of the field at a final time, 
defined as a multiple of tq . This approach to determin- 
ing the final kink number leads to slight discrepancies 
in the estimate of a in the extreme high/low-dissipation 
cases. These are nevertheless easily controlled 0, 
and as we will see below the accuracy obtained is more 
than sufficient for the needs of the present work. Finally, 
we perform a fit of the final defect density versus tq to 
a power law and thus obtain the scaling exponent a. 

In Fig. Q] we show the results for a case where the evo- 
lution before the transition is well in the under-damped 
regime. In 0] we showed that for this specific model, 
the transition between the two types of evolution takes 
place at about a c = 0.8. The two values of a< = 0.4, 0.5 
used here are both comfortably below a c . In both cases 



FIG. 2: Defect scaling power as a function of the noise am- 
plitude for t < 0. For t > the amplitude are a> = 0.4 and 
a> = 0.5. Error bars are as in Fig.0 



we allow a> to vary between 1.0 and 3.5, covering a part 
of the over-damped region of parameter space. As illus- 
trated in Fig. the results show clearly that the scaling 
parameter a is determined exclusively by the value of 
a > . The result for the scaling exponent is unmistakably 
typical of the over-damped regime, despite the fact the 
part of the evolution was under-damped. Not only that, 
but the precise value of a does not seem to affected by 
a < . Within the error bars, the curves obtained are indis- 
tinguishable from the case where the dissipation is kept 
constant, and equal to a>, throughout the transition. 

In Fig.|21we show the results for the opposite situation. 
Here we kept the dissipation after the transition constant 
and in the under-damped regime, while keeping a< in the 
over-damped region. The measured a is constant within 
error bars, and has a value typical of an under-damped 
evolution. The values of the scaling power for a quench 
with constant a set to 0.4 and 0.5 are a = 0.39 and 
a = 0.35 respectively . As in the previous example, the 
values shown in Fig. [21 coincide with these. Once again 
this strongly suggests that the outcome is determined 
exclusively by the period of the quench taking place after 
the phase transition. 

For our second example we take a more radical depar- 
ture from the original model @. While maintaining the 
dissipation at a constant value throughout the whole evo- 
lution, we change the rate of the quench at t = 0. For 
negative times the squared mass term m 2 (t) will go from 
/i 2 at t = —tq < to zero at t = 0. For t > an alter- 
native quench rate is introduced, with m 2 (t) decreasing 
linearly and reaching — ji 2 for t = tq>. As in the previ- 



1 Since we are working on the region of very low dissipations these 
are slightly larger than 1/3, the typical value for the under- 
damped regime, a consequence of saturation phenomena Il7tll8l . 
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FIG. 3: Defect scaling power as a function of the fixed, t < 
quench time tq< = 2 n< , for under-damped (a — 0.4) and 
over-damped (a = 2.5) cases. Error bars are as in Fig.0 



ous section, we expect the outcome to reflect the relative 
importance of the two periods of the evolution. 

In Fig. 03 we show the results for two cases, with a in 
the under-damped and over-damped region respectively. 
For both sets of simulations, tq < was kept constant at 
a value tq< = 2 n< for fixed n < whilst tq> = 2 n> was 
varied throughout the usual range, with n> = 1, ...9. 
If the final number of defects were determined by tq< 
we should observe no scaling, with a constant number of 
defects being obtained at the end of each simulation. In- 
stead, we found that the final density clearly scales with 
tq > , with a power that is determined by the magnitude of 
a. As seen in Fig.|3|the value of the scaling power changes 
very little with the quench rate at negative time. For the 
two choices of a — 0.4, 2.5, the exponents for the basic 
quenches with tq < — tq> are 0.39 and 0.23 respectively. 
The measured ct's lie very near these, confirming that the 
t > part of the evolution is the dominant one. As an 
extra check we simulated the reversed situation where n> 
was kept fixed and n< was allowed to vary in the usual 
range. As expected the final defect density was constant 
within error bars, showing no dependence on the value 
of tq< for negative times. Once again this highlights the 
fact that the mechanism of defect formation is dictated 
by what happens after the transition takes place, in vari- 



ance with some of the arguments of Q, 0, Q . 



V. CONCLUSIONS 

If there had been any doubt before, we have shown 
conclusively that defect density is controlled by the cool- 
ing of the system after the critical temperature has been 
passed. It is true that correlation lengths do freeze in be- 
fore the transition. However, their defining long-distance 
behaviour is irrelevant to defect densities, which are con- 
trolled by the short-distance behaviour that determines 
field zeroes. 

As we have stressed, transitions with discretely differ- 
ent behaviour before and after the critical temperature 
of the type that we have examined are unlikely to have 
direct physical counterparts. However, the transitions 
in Josephson Junctions provide a more subtle example 
where behaviour before and after the transition is very 
different. Before the conductor-superconductor transi- 
tion we have two bulk conductors, individually amenable 
to the Zurek analysis of freezing in of correlation lengths 
For these, is appropriate, with a < = 1/4. How- 
ever, after the transition we have the Josephson effect, for 
which field ordering is not determined by the behaviour 
of the individual superconductors, for which would 
have given an identical <t> = 1/4 3]. Instead, field or- 
dering is constrained by the velocity of light in the oxide, 
the Swihart velocity [22|]. Naively, we could not easily 
distinguish between the two possibilities since, for ide- 
alised symmetric junctions this again gives the same scal- 
ing behaviour a > = 1/4 j2^ after the transition, albeit 
with a different prefactor. However, for realistic junctions 
proximity effects and only partial critical slowing-down 
er< = 1/4 enforce cr> = 1/2 16]. Empirically, the data 
is described well by the latter [l(j , and is not compatible 
with cr< = 1/4. 

Both the analytical estimates and the numerical simu- 
lations presented above suggest strongly that our results 
should apply to domain formation in generic second-order 
phase-transitions. 
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